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ABSTRACT 

Outflows from active galactic nuclei may be produced by absorption of continuum ra- 
diation by UV resonance lines of abundant metal ions, as observed in broad absorption 
line quasars (BALQs). The radiation pressure exerted on the metal ions is coupled to 
the rest of the gas through Coulomb collisions of the metal ions. We calculate the 
photon density and gas density which allow decoupling of the metal ions from the rest 
of the gas. These conditions may lead to an outflow composed mostly of the metal 
ions. We derive a method to constrain the metals/H ratio of observed UV outflows, 
based on the Lya and Siiv AA1394, 1403 absorption profiles. We apply this method 
to an SDSS sample of BALQs to derive a handful of candidate outflows with a higher 
than solar metal/H ratio. This mechanism can produce ultra fast UV outflows, if a 
shield of the continuum source with a strong absorption edge is present. 

Key words: galaxies: active - quasars: absorption lines - quasars: general. 
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1 INTRODUCTION 

Broad absorption line quasars (BALQs) are a subclass 
of active galactic nuclei (AGNs) with spectra which ex- 
hibit broad, strong, and blue-shifted absorption features 
l|Reichard et al.l |2003| and references there in). Approxi- 
mately 15 per cent o f quasars are classi fi ed as BALQs 
llReichard et al.1 l2003l : iKnigge et al] l2008l : IScaringi et al] 
l2009t ). Although there are some differences in the contin- 
uum and emission lines between BALQs and non-BALQs, 
they appear to be drawn from the same pa rent population 
l|Wevmann et al.lll99ll : lReichard et al.ll2003l '). The broad ab- 
sorption lines (BALs) are the most prominent spectral fea- 
ture indicating an outflow of material driven from and by 
the AGN. This outflow is a compelling candidate for the 
feedback mechanism of the AGN on its host. The outflow 
is commonly invoked as a growth and activity regulator 
l|Di Matteo. Springel fc Hernauistll2005l and citations there- 
after) and as a mechanism of metal enrichment of the ISM 
and IGM (|Mo11 et al.ll2"007h . 

How do the outflows form? Although much theo- 
retical work has been done to model AGN outflows, a 
work which stems from studies of stellar winds (e.g., 
iLamers fc Cassinellil Il999l ). the driving physical mecha- 
nism of outflows remains uncertain. The first models 

involved only gas a nd radiation pressu re as a driving 

force o f the outflow llVitello fc S hlosmanl 19881: Aray fc Lil 
1 19941 ; lArav, Li fc BegelmanI (r994i : iMurrav et al.1 1 19951 : 
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iProga. Stone fc Drew! Il998l : IChelouche fc Netzeij 1200 ll ). 
Later wo rk introduced more detailed models of radiation 
pressu re (|Proga. Stone fc Drew! 19991: Chelouche fc Netzerl 
l2003al ) and magnetic flelds (|de Kool fc BegelmanI Il995l ). 
The effect of rad iation-drive n wi nd on the emerging 
spectrum (lArav fc Bege lm anI [l994l ) and galactic bulge 
l|Fabian. Celotti fc Erlundi 120061 ) has also been consid- 
ered. More recently, with the advance in compu- 
tational power, more detailed magnetohydrodynamical 
and magnetocentrifugal models of the outflow have 
been 



iigl fc Kartiel il993: jProgal 
Proga. Stone fc Kallma'nH200d: iProgal l2003l: lEverettl 



developed (iKonii 



200C; 



2005 



Kurosawa fc Progal 20081 ). Currently, the models pro- 
duce almost no predictions which can help discrimi- 
nate between them usi ng observations (see a review by 
ICrenshaw. Kraemer fc Geo rge 2003 ) . 

In current models the outflow is treated as a single- 
component fluid, and the radiation pressure is parametrized 
by a "force multiplier", which encapsulates the effects of 
radiation pressure due to Thomson scattering and line 
and dust absorption. What effect does the radiation pres- 
sure have on the differe nt constituents of the fluid? A s 
was first pointed out by ISpringmann fc PauldrachI (|l992l ). 
for radiation-driven stellar winds, a single-component fluid 
treatment overlooks the p ossibility of metal ions decoupling 
from the mostly-H outflow. ISpringmann fc PauldrachI (Il992l ) 
have shown that metal ions can run away (i.e., decouple) 
from the proton-electron fluid, if the coupling between ions 
and the fluid is small enough, compared to the photon flux 
absorbed by the ions. Photon absorption by a metal ion 
results in a net gain of momentum i.e., ion acceleration (as- 
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suming illumination by an anisotropic source) , while the ac- 
celeration of H is inefflcient as H is mostly ionized. The ion is 
decelerated by Coulomb scatterings by the ambient protons 
and electrons. The deceleration is velocity dependent and is 
maximal when the ion has a velocity similar to the thermal 
velocity of each of the constituents of the proton-electron 
fluid. The ion can run away from the fluid, if the rate of pho- 
ton absorption (i.e., flux density) is large enough, so that the 
accumulated gain in velocity of the ion surpasses the fluid 
thermal velocity before the ion is significantly deflected by 
Coulomb scatterings i.e., if the acceleration due to the ra- 
diation pressure is larger than the deceleration due to the 
Coulomb force. The metal decoupling prevents the proton- 
electron fluid from accelerating further and the decoupled 
metals will produce a fast wind through the proton-electron 
fluid. 

Metal ion decoupling from the mostly-H gas may be 
relevant in AGN outflows due to two reasons. First, gas in 
close vicinity of the active nucleus is exposed to a large flux 
of radiation; second, strong ion absorption is observed. Our 
goal is to explore whether the occurrence of metal decoupling 
is a plausible process in AGN. Self-consistent models of an 
outflow and the structure of the ionized gas region, which 
incorporate the process of metal ion decoupling, are beyond 
the scope of this study. 

How can an occurrence of the metal decoupling be ob- 
servationally verified? The metal decoupling will mostly af- 
fect the absolute metallicity of a wind i.e., metal to H abun- 
dance ratio. The decoupling effect on the relative i.e., metal 
to metal abundance ratio, might be negligible. Since the 
wind can be composed from several types of decoupled met- 
als, the relative abundance ratio of the wind might be similar 
to that of the wind origin, where the metals are mixed with 
H. Thus, indirect (i.e., relative) outflow metallicity measure- 
ments cannot be used to test the proposed metal enrichment 
scenario. Direct (i.e., absolute metals to H) metallicity mea- 
surements are needed to test the validity of the enrichment 
scenario in AGN. 

There is accumulating evidence of metal abundances 
highe r than solar for the broad line region (BLR; iHamannl 
1 19971 . iHamann fc Ferlandl Il999l ). There is also evidence 
for higher than solar metal abundances in AGN outflows 
from analysis of narrow absorption lines (jArav et al ] |2007l : 
IWu et al.ll2O10l: IHamann et al.|[2oTll ) and broad absorption 
lines ( Hamann et al.ll 19971). However, fo r some objects the 
results are inconclusive (jHamann 1998) or consistent with 
solar abundances (|Arav et al. 20011. The main difficulty in 
estimating the element abundance arises from a large uncer- 
tainty in the measurement of ionic column densities when 
the covering fact or is velocity depend ent (iHamannl Il99^ : 
lArav et al.lll999l ). Note however that ICottis et air(j2010l ) 



disfavour radiation-driven outflow altogether, due to the ab- 
sence of an excess of objects displaying line-locking between 
Lya and Nv. A metallicity of a few times solar can be ex- 
plained b y stellar enrichment due to extended star forma- 
tion (e.g.. lHamann et al.ll2002|). as sug gested in high-redshift 
quasars (z > 3.5. iDietrich et al.ll2003i ). However, if a metal 
ions runaway occurs, the metal/H column ratio can be larger 
by orders of magnitude, compared to the solar abundance 
ratio. 

Below we derive the conditions for an ion runaway in 
gas in the vicinity of an active nucleus. We also present a 



method to set a robust and direct lower limit on the metal/H 
ratio in outflows based on UV absorption lines. The method 
is then implemented on a sample of AGN f rom the Sloan 
Dig ital Sky Survey (SP SS; lYork etahlbOOOl ) BALQs cata- 
log (|Scaringi et al. |2009| ) . The ion runaway model is outlined 
in Sec. 2. The direct metallicity estimation method and its 
implementation are described in Sees. 3 and 4, respectively. 
We discuss the results in Sec. 5. Our conclusions are sum- 
marized in Sec. 6. 



2 METAL IONS RUNAWAY 

We begin by describing the conditions for the separation of 
metal ion from an ionized H gas (i.e., runaway). First, we 
estimate the photon absorption rate by the metal ion, and 
the mean momentum gain by the absorption of each photon. 
A product of the two quantities is equal to the accelerating 
force on the ion due to radiation pressure. Second, this radia- 
tive accelerating force is compared to the decelerating force 
that is caused by Coulomb scattering of the metal ion by the 
H gas (i.e., protons and electrons). Finally, we investigate 
the conditions (photon flux density, H number density and 
temperature) that lead to an accelerating force larger than 
the decelerating one. We first present a qualitative order- 
of-magnitude estimates, and then a more exact and formal 
derivation of the separation mechanism. 



2.1 Qualitative estimates 

An estimate of the ionization parameter \J for a metal ion 
runaway in AGN can be roughly estimated as follows. The 
mean momentum gain of an ion due to absorption of a pho- 
ton is 



h 



(1) 



where h is Plank's constant and Aion is the wavelength of 
the absorption transition. The mean increment of the ion 
velocity due to photon absorption is Itv = p/mion, where 
mion is the ion mass. A typical value for an UV absorption 



line, e.g. C iv, is 
rate is 



20 cm s . The photon absorption 



r^i = BuJ, (2) 

where B12 is the Einstein absorption coefficienlQ, and J is 
the mean intensity averaged over the metal ion absorp- 
tion line profile i.e., J = Jv4>(i')dv, where </'(i^) is the 
line profile function (e.g., [llvbicki fc LightmanI lioOll ) . The 
value of B\2 for typical resonance UV lines (e.g., Civ) is 
~ 7 X 10^ cm'^ erg~^ s~^. We begin by estimating r-^ in the 
BLR in AGN. We assume a bolometric luminosity of Lboi = 
10^^ erg s"^ and vL,, = 2.5 x 10'*'^ erg s~^ at A = 1500 A, 
which yields L^(1500 A) = 1.3 x 10^° erg s"^ Hz'V The 
implied BLR distance from the central illuminating source 
is i?BLR = 0.1 pc, and the resulting mean intensity is J ~ 



^ The Einstein absorption coefficient is given by iJi2 = 
iTT^e^\a^fl2/hmeC^ , where e and me are the electron charge 
and mass, /12 is the transition oscillator strength and c is the 
speed of light. We adopt the values of /12 and Ajon from lMortonI 
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J„(1500 A) ^ 10"^ erg s"^ cm-^ ster"^ Hz~\ Thus, the 
photon absorption rate is ~ 7 x 10^ ■ 1 x 10"^ « 10^ s"^. 

The decelerating force acting on the metal ion is due to 
Coulomb collisions with the ionized H gas constituents i.e., 
protons and electrons. A simplified condition for a runaway 
is that the metal ion gains a velocity larger than the thermal 
velocity of the ions «th,ion between consecutive collisions. 
Mathematically this condition can be stated as 



(3) 



where rcoi is the collision rate. Noting that the main contri- 
bution to the decelerating force at low ion velocities origi- 
nates from collisions with protons (see Sec. l2.2|) . we consider 
collisions only with protons for the remainder of this quali- 
tative description. The collision rate is 



^col = WpHth.pCTcol, 



(4) 



where rip is the number density of the protons, Vth,p is the 
proton thermal velocity and Ucoi is the cross-section for a 
proton-ion collision. The collision cross-section can be esti- 
mated by evaluating the impact parameter b at which the 
metal ion is scattered by 90 deg. The 90 deg scattering oc- 
curs when the Coulomb potential energy and the kinetic en- 
ergy of the ion (Skin) are equal i.e., Zionfi^/b ~ i?kin, where 
e is the electron charge and Zion is the ionization state of 
the ion (~ 3). We take i?kin to be the thermal energy i.e.. 



kT 



10" 



erg for a typical BLR temperature 



of r = 10 K, where k is the Boltzmann constant. This 
yields fe ^ 5 x 10"^ cm and ctcoI = Trfc^ ^ 8 x 10~^^ cm^. 
The ion and protons thermal velocities for T — 10* K are 
Vth.ion « 4 X 10^ cm s~^ and «th,p « 1 x 10® cm s~^; and 
from Eq.El rcoi < 20 x 10^4 x lO'^ ~ 5 x 10"^ s"^ Now we 
can find the condition on the number density from Eq. 21 
Wp < 5 X 10"V(1 X 10® ■ 8 X 10"") ~ lO'^ cm"^ in order to 
get a runaway acceleration. 

Adopting the accelerating and decelerating forces es- 
timated above, one can evaluate Ur = n-y/rip that pro- 
duces a runaway, where is the ionizing flux density. Since 
rip = 10^" cm~^ is the typical [/ = 0.1 BLR density, one 
needs = 0.1 x (lO^'^/lO^) 10* to have a runaway. The 
value of Ui is independent of R, since r-^ oc , and thus 
rip oc , and also n-, cx R~^ . At this large value of Ui, 
photoionized gas is fully ionized, and the relevant UV ab- 
sorbing metal ions are not present. A runaway thus requires 
a shield between the illuminating source and the gas. The 
shield has to filter the high energy photons that can overi- 
onize the metals (e.g., > 4.7 Ryd for C-outflow) Note that 
such a shield appears to be required to explain radiation- 
driven BALQ-outflows fe.g.- lMurrav et al.l[l995l l. 



2.2 Quantitative estimates 

We describe a microscopic model in which an ion (e.g., C'^^) 
is embedded in a fully ionized gas mostly made of ionized 
H, illuminated by a continuum in the radial direction. In 
addition, we make the following simplifying assumptions: 

(i) The ion is assumed to have only one transition i.e., 
only one absorption line at the laboratory wavelength Aion. 



(ii) All calculations are conducted in the non-relativistic 
limit. 

(iii) The gas is in a free-fall circular motion, held by grav- 
ity and radiation pressure. 

The mean momentum transformed from a photon to an ion 
in a process of absorption followed by reemission is given 
by Eq. [T] The photon absorption rate (Eq. [2} expressed in 
terms of oscillator-strength /12 and flux density (note 
that Ji, = fv/An) is 



7re 



-/l2/i/(Aio 



T — " 1 

meC/lc/Aion 

The acceleration due to radiation pressure is 



arad = P/ m-ion ■ = 



T^e^ /12 



/.(Aio 



(5) 



(6) 



The main relevant UV transitions, C iv, Si iv and N v, have 



similar atomic parameters (/12/mion 
10"74Aion) [cm s-2]. 



10^^ gr~7 yielding 



flrad 



(7) 



We adopt the met hod described by 
ISpringmann fc PauldrachI (|l992l ) for multicomponent 
fluid winds driven by radiation of hot starslf] Protons and 
electrons induce a deceleration force on the test-particle 
due to Coulomb friction (e.g., ISpitzeJ Il962h . The total 
deceleration is 



flfric = w» ^,^ lnAG(a:i 



Z — p,G 



fcTmi, 



(8) 



where rii is the number density of constituent i (note that 
rie ~ rip ~ wh, where the latter is the H number density), 
Zion is the ionization state of the test particle. In A is the 
Coulomb logarithm, G{y) = [^(y)— 2/$'(y)]/2y^ is the Chan- 
drasekhar function with ^{y) being the error function and 



mionmi 



JTlion + rUi 



(9) 



where w is the relative veloc ity between the test particle and 
the fluid fe.g.. lSpitzerlll962l). For the C oulomb logarithm we 
adopt the equation from ISpitzerl (|l962l ). 



InA = In- 



Trnn 



(10) 



1 Ryd = 13.6 eV 



Several authors use the total number den sity of free particles 
instead of riH (e.g.. lOwocki fc Pulj|200"3 ). This has a minor 
effect on InA, lowering it by few per cents. 

Figure [1] presents the Coulomb friction deceleration 
as a function of the relative velocity 10, while assuming 
wh = 10® cm"^, T = 10* K and a C^+-like test-particle 
ion. We plot both the individual decelerations due to pro- 
tons and electrons, and the total friction deceleration. Pro- 
tons contribute to Ofric mostly for w < 10^ cm s^^, while 
electrons contribute above those velocities. 

The method can be derived fro m the "classical" description 
of Coulomb friction llSpitzeJll96l) by adopting a reduced mass 
instead of the field mass in the equation for velocity normaUzation 
(denoted by Zf in Spitzer). This adaptation is needed because the 
derivation in Spitzer is made for an electron and not for an ion 
particle. The reduced mass with a proton is « rrie for the former, 
but it is nip for the latter, where nip is the proton mass. 



© 2012 RAS, MNRAS 000.[TlfT6l 



4 A. Baskin and A. Laor 



The maximal frictional deceleration [a'-l^^^) for each 
constituent is reached when w ~ «th, where «th is the 
thermal velocity of the constituent. This is because G(xi) 
peaks aX xi ^ I [G(l) = 0.214 (e.g.. iKulsrudI 120051 ')]. and 
JTiion ^ me,mp which yields Xi « w/y^2kT/mi = w/vth,i- 
Once a test particle reaches w > vth.c, the total frictional 
decelerating force on it diminishes and it runs away (if the 
accelerating force remains active). Since the dependence 
of In A on T and nu is logarithmic (Eq. llOp . and T and 
riH ^ 1 for an AGN environment, InA « const. This ap- 
proximation simplifies the a™/^" dependence on T and nn 
to a™^" oc uhT'^. The a'^^^ can be estimated to a good 
approximation by 

loga'St^ ^ {S.Sti.D + log ^ [cm s'^], (11) 

where the range on the constant is the maximal variation 
of Eq. m [taking G = G(l)] for the C^+, N*+ and Si=^+ ions 
and In A, when riH is varied in the 10^-10^^ cm~"^ range 
and T in the 10^-10*' K range. Note that for T = lO'^ K 
Vth,e ~ 0.02c and raising T to much higher values will break 
the non-relativistic approximationQ 

An ion will run away when a^d > cif^it, i.e. when 

log/4Aion)-lognH + logT> -6.2. (12) 

Assuming a spectral energy distribution (SED), we can con- 
vert Eq. [T2]to a condition on U. An SED with fi, oc 
up to 1 Ryd, and a steeper slope beyond that, gives n-, ~ 
vf^/hvc Si 5 X 10^V-^(1 Ryd) « 3 X 10^V>'(Aion). Inserting 
in Eq. [12] yields 

log - 15.5 - log nn -HlogT > -6.2 (13) 

or simply, log 17 > 5.3 for logT = 4. 

The model described above and the derived conditions 
are only estimates of the conditions needed to produce a 
runaway of metals. One needs to sum up over all resonance 
lines. Also, the ionization state of the metal ions changes 
with time, and the relevant resonance lines change. For ex- 
ample, a C atom is found ~ 60 per cent of the time in the 
C^^ ionization state, when the latter is the dominant state 
(e.g.. Fig. [21 top panel). Note that the electric field produced 
by the separation of a metal ion causes the ion to "carry" 
with it the proper number of electrons, and the gas remains 
globally neutral. 

The difference in ionization fractions of metals com- 
pared to H is the main reason why the radiative acceler- 
ation is more efficient for metals than for H. When the gas 
is ionized so that there is a significant fraction of elements 
that are ionized three and above times (C^^, Si^^, N'*''', 
etc.), H atoms ar e mostly ioniz ed and the fraction of H" is 
- IQ'^U'^ (e.g.. lFerlan"dlll999l : see also Fig. [2] top panel). 
The for H'' and metal ions is comparable, but the metal 
ions spend most of the time at the dominant ionization state, 
while H spends only a fraction of 1Q~^U~^ of the time as 

Note that for a very large temperature (T > 10® K) one needs 
to use the de BrogUe wavelength Ajb a s the distance of minimal 
approach in the calculation of A fe.e.. [KulsrudlfioOSh . Even for 
T > 10^ K one needs to use Ajb for interactions with electrons. 
This introduces a small negative correction to In A, which causes 
'^fric^e ^ ''fric'p' Thus, disregarding this effect does not affect the 
results of our study. 
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Figure 1. The Coulomb friction deceleration as a function of 
relative velocity between a test-particle and a fluid. The decelera- 
tion due to protons (dashed line), electrons (dotted line), and the 
total deceleration (solid line) are presented for = 10® cm"'', 
T = 10"' K and a C^+-like test particle (Eq. Ell. The radiative 
acceleration of a C''"'" particle absorbing at Civ A1549 transition 
(dot-dashed line) is also plotted for two different /y(1549 A). Note 
that the maximal deceleration o f each fluid constituent is achieved 
when w K, iith, where Uth = \/ 2kT jm-p^^ is the thermal velocity 
of the constituent. An ion will run away when a^ad > ifric- 



H , and being accelerated. Therefore, the relative flrad of H 
is reduced by a factor of 10~^C/~^ and is negligible. 



3 ESTIMATING H/METAL COLUMN RATIO 
IN OUTFLOWS 

The straightforward method to directly estimate the abso- 
lute metallicity of an outflow is to measure the ionic column 
density ratio of metals to H. The difficulty to implement 
this method for quasars is that usually the H column den- 
sity, derived based on Ly a absorption, is poorly constrained. 
Most spectra are obtained from ground based observation of 
> 3 quasars, which are subject to significant intervening 
Ly a absorption systems. This foreground absorption makes 
it difficult to determine the unabsorbed emission level. In 
addition, the potentially absorbed blue wing of Ly ol cannot 
be estimated based on the unabsorbed red wing, as the red 
wing is blended with the Nv emission line and is also of- 
ten partly absorbed by the Nv absorber. Thus, we choose 
to constrain only the minimal allowed A''(H'') in a statisti- 
cal sense by utilizing composite spectra of non-BALQs (see 
Sec.[4ll. 

The structure of this section is as follows. First, we 
search for the best metal absorption line that can be used 
to place a direct lower limit on the gas absolute metallicity 
(Sec. I3.l"|l . Then, we adopt photoionization models to set a 
lower limit on the necessary Lya absorption from that gas 
(Sec. 13. 2p . Finally, we present the procedure for estimating 
a lower limit on the absolute metallicity i.e., the H/metal 
abundance ratio in the outflowing gas (Sec. I3.3]l . 
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3.1 Which absorption line to use? 

To directly constrain the H abundance associated with a 
metal absorber, we need a UV metal line which is expected 
to be associated with a large A'^(H'') column. Three criteria 
are imposed to find the best metal absorption line. 

(i) The line can be observed with Ly a at the same wave- 
length coverage of a typical instrument (e.g., the SDSS cov- 
ers the 3800-9200 A range, and requires z > 2.7 to cover 
the Lya-Civ A1549 wavelength region). This criterion dis- 
favors Mg 11 A2799 which is found too far to the red relative 
to Lyof. 

(ii) The line should be prominent and not blended with 
absorption lines of other metals. This criterion disfavors the 
remaining low ionization UV lines. Although C+ (C ll A1335) 
and Al"*"^ (Aim A1857) produce larger column density ratios 
N {SP) / N {ion) than the higher ionization ions, their absorp- 
tion profile is generally shallower, and the implied minimal 
Lya absorption can be harder to detect. These lines are 
therefore excluded. 

(iii) The photoionization models should imply the largest 
minimal Lya optical depth r for a given r of metal line. 
Since r oc /i2Aion x A'^(ion), this is equivalent to requiring 
a minimal A'^(H°)/A''(ion) x l//i2Aion as high as possible. 
The value of Ai'(H'')/A''(ion) is calculated by photoioniza- 
tion models, and /i2Aion is set by the atomic physics. This 
criterion favours Siiv (see Sec. I3.2.2|) . which is observed to 
be as prominent as C iv in absorption for the objects studied 
here. 

Thus, the above criteria favour by elimination the Silv ab- 
sorption line. 



3.2 The photoionization models 

What is the minimal possible A'^(H'') for a given A'^(ion)? 
We estimate this quantity using the photoionization code 
CLOUDY for a very wide range of model parameters. Cal- 
culations are perfor med with version 08.01 of CLOUDY 
(jPerland et al.l[r998l ). A plane parallel geometry is assumed. 
The gas is assumed to have a solar composition. The ion- 
ization parameter U and the total H column density E are 
varied by 0.5 dex in the following range: —3 ^ log [7^2 
and 15 ^ log E ^ 24. The H density rip is varied by 1 dex 
in the range 5 ^ log Wp ^ 12. It should be noted that 
the assumed Wp has a small effect on the deduced minimal 
iV(H")/iV(ion). Models with 3 ^ logC/ sC 5, logE = 25-26 
and log rip = 8 are also calculated to find an absolute mini- 
mum on A'^(H°)/A^(ion) (see below). Note however that for 
these parameters the gas is optically thick to electron scat- 
tering (re. sc. > 1), and CLOUDY results may be incorrect, 
since the code is not designed to simulate regimes with Te.sc. 
well above 1. 



3.2.1 The spectral energy distribution 

The minimal possible value for A^(H'')/Ai'(ion) is mostly set 
by the assumed SED of the illuminating source. The value of 
A''(H°)/A^(ion) at a given U and rip is set by the number of 
photons able to ionize H" relative to number photons which 
produce the ion without ionizing it to the next ionization 
state. Since H° is mostly ionized by photons with energies 



at ~1 - 2 Ryd, a steeper SED in the 1-5 Ryd range implies 
more H'' ionizing photons for a given ionic column of Si^^ 
and C*"^ (for a fixed rip). Thus, a steeper slope produces a 
lower iV(H°)/Ar(ion). 

We adopt a broken power-law SED for the incident con- 
tinuum. We use the steepest SED consistent with observa- 
tions to produce the minimal A'^(H")/7V(ion). The adopted 
spectral slopes a (/^ oc u'") are described below. In the 
range 3000 - 10 ,000 A, we use the slope measured by 
lElvis et all \imi ) (a cutoff is assumed above 1 /im). In the 
UV range (1000 - 3000 A), we use the slope which is mea- 
sured using the flux of the composite BALQ spectrum (see 
Sec.gT]) between ~ 1300 and 2050 A. The measured a = 1.8 
is extremely red compared to the typical AGN a — 0.5 in 
this wavelength range, but since the N(R°)/NiSi^+) de- 
pends on energies larger than 1 Ryd (i.e., A < 912 A), the 
results of this study are not sensitive to the precise value 
of a below 1 Ryd. In the 500 - 1200 A range, iTelfer et al] 
(2002) find a of 2.0 and 1.6 for radio-loud (RL) and radio- 
quite (R Q) quasars, respectively. Although m ost BALQs are 
RQ ()Stocke et al.lll992l : iBecker et al]|2000l 'l. we adopt the 
steeper slope of RL quasars for A = 500—1000 A, because we 
study carefully selected objects which may happen to be RL. 
There are no direct meas urements of a in th e ~ 60 — 500 A 
range, and we extend thel^fereFaD (|2002l ') slope down to 
12 A (1 ke V). Note that a = 2 is steeper than the slope 
reported bv lLaor et al.l (| 19971 ') for the soft X-ray range (0.2 
- 1 keV; a = 1.2 and 1.7 for RL and RQ quasars, respec- 
tively), but this energy range has lit tle effect on the results of 
our st udy. In the 2-10 keV range. iBrandt. Mathur fc Elvid 
(|l997l ) measure a = 1, and iReeves fc Turneil (|200G| ) report 
a = 0.7 and 0.9 for RL and RQ quasars, respectively. We 
adopt Q = 1 for the whole hard X-ray range (1 - 100 kcV), 
and assume a cutoff above 100 keV ( Molina et al.. ,20091). 
The overall adopted SED has Qox = 1.9 between optical 
(3000 A) and X -ray (2 keV), whi ch is close to the typical 
BALQ Qfox = 2 (jBrandt. Laor fc W ills 2000). Table [1] sum- 
marizes the adopted slopes and the respective energy and 
wavelength ranges. 



3.2.2 The model results 

Figure [5] presents the ionization structure of the slab (top 
panel), and the integrated A'^(H°)/A'^(ion) as a function 
of integrated A''(ion) (bottom panel; this is equivalent to 
assuming different E), assuming log 17 = —1 and rip — 
10** cm"''. Although the Si^"*' ionization fraction is larger 
by ~ 3.5 dex than H" in the region where the Si^^ has a 
significant fraction and column (top panel), the integrated 
iV(H")/7V(Si^+) ~ 10 (bottom panel) because H is ~ lO^ '^- 
time more abundant than Si (assuming Zq). Note that C^^ 
and N'*"^ place a lower constraint on A'^(H'') than Si^"^. A 
higher constraint is placed by C^. But, the observed ab- 
sorption profile of Cll A1335 is shallower than Silv A1397 
(see Sec. |4]) and yields a lower overall constraint on Ly a ab- 
sorption. Thus, we chose Si^^ as the optimal ion to constrain 
Af(H°). 

Figure |3] presents the dependence of iV(H°)/iV(Si^+) on 
Ai'(Si^^). The dependence is plotted for different values of U 
assuming rip = 10* cm~^. For log (7 = —2 we also plot the 
dependence for rip = 10^ and 10^^ cm~^, which produces 
similar curves. The important result of the models is that 



© 2012 RAS, MNRAS 000.[TlfT6l 



6 A. Baskin and A. Laor 



Table 1. The adopted spectral energy distributional 
A rangJS 

(A) 



Energy range 
(eV) 



o|Sl References 



1.2 - 4.1 
4.1 - 12.4 
12.4 - 24.8 
24.8 - 1000 
(1 - 100) X 10^ 



3000 - 10,000 
1000 - 3000 
500 - 1000 
12 - 500 
0.1 - 12 



"05 lElvis et alJ lll994ll 
1.8 this papeiBl 
2.0 iTelfer et al] | |2002|) 

2.0 ^ 

1.0 Brandt et al] \l99t ): 
, Reeves ^"^irneil 
1I2OOOIUI 



" A cutoff is assumed below 1 and above 100 keV. The overall 

SED has a slope of Oox = 1-9 between 3000 A and 2 keV. 

^ The corresponding wavelength range for the listed energy 

range. 

oc u-'^. 

The slope is the steepest one consistent with the composite 
BALQ spectrum between ~ 1300 and 2050 A (see text). This 
slope is then extrapolated to 1000 and 3000 A. 
^ Extrapolation of the Telfcr ct al. (20o|) slope. 
/ Note that ISrandt et al.l li 19971) and iReeves fc Turned 1I2OOOI) 
measurements were conducted in the 2-10 keV range. 



N{'H°)/N{Si^+) > 10 for 10^^ < N{Si^+) < W^^ cm-^; for 
lower values of A'^(Si^''') the minimal value of the ratio goes 
up to ~ 50. When A''(Si^+) > 10^*^ cm"^ the ratio drops a bit 
below 10, but for that range iV(H'')> lO" cm"^ and yields 
a saturated Ly a absorption profile. Raising U above 10 does 
not lower the minimal Af(H°)/Af(Si iv) below 0.4 dex. This 
absolute minimum of logiV(H°)/iV(Siiv) = 0.4 is achieved 
even for an extreme high value of log (7 = 5.0 

Using the derived values of A''(H'')/A^(Si'^^), one can set 
a robust lower limit on r of Ly a based on the measured r 
of Si IV. The ratio between the two optical depths is set by 
the ratio of /12A of the two lines and by N(}i°)/N{Si^+). 
This yields the robust lower limit of 



r(Lya) ^ 0.71 



iV(HO) 



iV(Si3+) 



X r(Silv A1394), 



(14) 



wher e the numerica l value is calculated using /12A values 
from iMortoril (|l99ll ). Note that [iV(H°)/iV(Si^+)]min is a 
function of Af(Si^+). The Lya optical depth is ~ 7x r(Silv) 
for the typical value of [Ar(H'')/Ar(Si^+)]min ~ 10. This 
means that Lya absorption is deeper than Siiv, unless 
T ^ 1 and the saturated absorption profile is set by a 
velocity-dependent covering factor f or bo th lines, and is 
therefore identical. Note that iKwanI l|l990l ) reached a sim- 
ilar conclusion that T(LyQ) > T(Silv) even while using a 
steeper slope in the EUV than used here; and while filling 
the Lya trough by scattered and emitted radiation from 
other regions. 



^ Note that log U = 5 requires log S 
log A''(Silv) > 17. This U and E imply Tc.sc. 
results may be incorrect. 



= 26 to achieve 
10, and CLOUDY 




10.5 11 11.5 12 12.5 13 13.5 14 

Log depth [cm] 



14.5 



.2 

bO 
O 




14 15 16 

Log W(ion) [cm^^] 

Figure 2. The slab ionization structure (top) and the ratio of 
column of H'' to the column of different ions versus the ion col- 
umn (bottom). We assume logC/ = —1 and rXp = 10* cm~^. The 
assumed abundances relative to H of C, N and Si are —3.61, —4.07 
and —4.46 in log-scale (CLOUDY default solar abundances). Al- 
though C+ can place a higher constraint on the allowed Af(H'^) 
than Si''"'" , it is not used in this study, because the observed C n 
absorption is weak. 



3.3 The procedure for setting a robust lower limit 
on Af(H") and the absolute metallicity 

The above described lower limit on r(Lya) for a given 
r(Silv) allows to place a robust lower limit on the absorber 
metallicity according to the following procedure. 

(i) Derive the intrinsic r(Silv) and Af(Si''"'') from the ob- 
served SilvAA1394, 1403 absorption profile. The derivation 
of the intrinsic r(Silv) depends on the assumed physical 
properties of the absorber as further described in Sec. 4.3. 

(ii) Find the minimal allowed iV(H°)/iV(Si^+) by pho- 
toionization models and derive the minimal A''(H''). 

(iii) Calculate the minimal Ly a absorption profile by the 
absorption system from step (i) using Eq. 1141 

(iv) Reconstruct the unabsorbed spectrum using the min- 
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14 
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Log iV(Si3+) 
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Figure 3. The column A''(H°) relative to the column A''(Si^+) 
versus the latter for different ionization models. The assumed 
ionization parameters are indicated. All models are calculated 
with rip = 10* cm~^ (solid lines). Two additional models with 
rtp = 10^ and 10^^ cm~'^ are plotted for log (7 = —2 (dot-dashed 
and dotted lines, respectively). Note the similarity between the 
three curves. The range of measured A'^(Si^+) in this study for 
a single absorption system is marked by two thick vertical bars. 
Note that in this range Af(IfO)/Af(Si^+) > 8 for all possible mod- 
els. 



imal Lya absorption i.e., divide the spectrum by the profile 
from the previous step. 

(v) If the reconstructed spectrum is consistent with a 
"normal" AGN spectrum, then the gas metallicity is consis- 
tent with solar. If an artificial emission feature is produced, 
continue to the next step. 

(vi) Scale down the minimal allowed A''(H'') from step (ii) 
by a parameter s, and repeat steps (iii)-(iv). Iterate until 
an "acceptable" AGN spectrum is achieved. The lower limit 
on the absorber metallicity is then s x Zq. 

This procedure can be trivially extended to cases where 
several prominent absorption systems are needed to recon- 
struct the Si IV absorption system. Note that the last step 
assumes that the minimal allowed N{}l'') / N{Si'^^) scales lin- 
early with the metal abundance. The assumption certainly 
breaks for metallicity larger than a few times solar due to the 
non linear response of the slab ionization structure with in- 
creasing metallicity. This does not hinder the identification 
of objects with super solar metallicity {Z > IOZq). But, 
placing a correct lower limit on Z for these objects requires 
a grid of photoionization models with Z > WZq. 



4.1 The BALQ data set 

The dat a set is drawn from t he SPSS DR5 BALQ cata- 
logue of IScaringi et al.l (|2009t ). IScaringi et al.l ()2009l ) con- 
structed the most complete catalogue of BALQs in DR5 (in 
the 1.7 < z < 4.2 rage). They cross-ch ecked the object clas - 
sification by their method with that of I Gibson et al.l l|2009l ). 
and visually inspected and classified the objects for which 
the two classifications disagree. This procedure lowered the 
probability for a false nega,tive c lassification from ~ 20 per 
cent for the Gibson et al.l ||2009|) catalogu e to ^ 4 per cent 
(|Scaringi et al.ll2009l ). The IScaringi et al.l l|2009h catalogue 
comprises 3552 objects. The object spectra were retrieved 
from the SDSS DR7. A successive series of selection criteria 
is imposed on the spectra. Each criterion is d escribed below 
in the same order it is imposed. Note that the lScaringi et all 
(2009) catalogue includes objects with BAL that are not 
shifted from the emission-line center. Since we want to min- 
imize the uncertainty on the measured absorption profile 
due to the underlying unknown emission profile, we require 
a blueshift of at least ~ 1000 km s~^ of the absorption from 
the emission-line center. The second criterion described be- 
low happens to exclude objects that do not meet this re- 
quirement. 

(i) We include only objects with z > 2.7 to ensure that 
the spectrum contains a sufficiently large (~ 200 A) wave- 
length region blue- ward of the Ly a emission line. The re- 
gion is chosen such that it is not affected by the intrinsic 
Lya absorption. This region is used in the following steps 
to constrain the intrinsic unabsorbed spectrum. This crite- 
rion leaves 973 objects in the data setlj 

(ii) Since the spectra have relatively low S/N (~ 5 per res- 
olution element at continuum regions unaffected by absorp- 
tionfl the spectra are smoothed by a 22 pixel-wide moving- 
average filter to achieve S/N of at least ~ 15 per bin [the 
width is equivalent to ~ 9 resolution elements (|York et al.l 
I2OOO!)]. We then inspect the spectrum by eye for objects hav- 
ing a Si IV absorption trough with width larger than 15 A 
(~ 3000 km s~^) but smaller than ~ 40 A, to prevent a sig- 
nificant blending between Ly a and N v absorption troughs, 
and normalized depth smaller than 0.5, to select objects with 
a prominent Siiv absorption. This yields a data set of 139 
objects. 

(iii) The Ly a and Si iv wavelength regions are examined 
by eye to evaluate the S/N within the absorption trough of 
each spectrum. Objects with poor S/N (< 1 per bin) in one 
or both regions are excluded from the analysis, leaving 78 
objects. 

(iv) In this step, we assume that the whole absorption 
profile is set by a velocity dependent covering factor (CP; 
i.e., a saturated profile). Note that the assumption yields the 
minimal possible absorption of Lya (see Sec. [S]). The Silv 
absorption profile is measured. This profile is used to re- 
construct the unabsorbed spectrum in the Lya wavelength 



4 SDSS BALQ SAMPLE ANALYSIS 

We investigate below whether there are any observed 
BALQs in the SDSS database that present evidence for a 
high metal/H outfiow. Derivation of the BALQ data set, 
construction of a non-BALQ control sample, spectral anal- 
ysis and the metal/H estimate are described below. 



° The redshift of 13 DR5 objects has been changed to values 
lower than 2.7 in DR7. Note t hat while the D R5 rcdshifts are 
used to choose objects from the lScaringi et al] (2009) catalogue, 
the revised DR7 redshifts are adopted to calculate the rest-frame 
wavelengths throughout the analysis. 

^ The S/N is estimated by dividing the mean fx by the standard 
deviation of fx in the Arcst = 1600 — 1640 A window. 
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region. Then, the following procedure is carried out to de- 
termine whether the reconstructed spectrum in the Ly a re- 
gion can he consistent with an unabsorbed Lya emission. 
A global power-law continuum is fitted for each spectrum 
using the measured f\ at ~ 1300 and 1700 A, where f\ is 
the observed flux density|f| Then, a synthetic Ly a emission 
line is placed on top of the continuum fit as an estimate of 
the intrinsic unabsorbed Ly a emission. Since the goal of this 
analysis is to find objects that require metallicity larger than 
solar with a high statistical significance, we use a synthetic 
line that has a large EW of 120 A and its FWHM is varied 
between 2000 and 10,000 km s"^ The chosen EW corre- 
sponds to the upper ~ 15 per cent of the Lya-|-Nv blend 
EW d istribution values, derived by lOiamond-Stanic et al.l 
l|2009[) for « > 3 non-BALQs. Objects for which the recon- 
structed spectrum is well above the synthetic emission lines 
(as determined by eye) are selected. This step leaves 13 ob- 
jects in the data set. 

(v) Spectra of 11 of the objects are similar, having a 
prominent presence of low-ionization absorption lines (C ll 
and Aim) in particular. We chose to analyse only those 
objects, since they might represent a particular subgroup of 
BALQ. The other two objects have a Si iv absorption trough 
shallower than Civ, and a weak Aim absorption. 

(vi) To make the sample more homogeneous in terms of z 
and Mi values, we exclude 3 more objects. The homogeneity 
is important in order to allow a comparison to non-B ALQs 
(see below). We use M, from ISchneider et all (|2oTol as a 
proxy for the luminosity, where i is the SDSS i filter|j Two 
of the objects have a relatively high z (~ 3.8) and one ob- 
ject has a high absolute magnitude (Mi ~ —28.7) and are 
therefore excluded for the sake of homogeneity. The final 
data set is comprised of 8 objects with 2.70 < z < 3.07 and 
-27.72 <M,< -26.18. 

Table[2]lists the name, z and h'U of the final data set objects. 
Note that no conclusion with a statistical meaning is drawn 
from the few objects in the final data set. The goal of the 
analysis is to check whether there are objects consistent with 
a high metal/H outflow. 

Figure |4] presents a composite spectrum of the 8 objects 
of the data set. Each spectrum is smoothed as described 
above [step (ii)] and normalized using a mean fx in the 
wavelength window of ~ 1600—1640 A. The composite spec- 
trum is a geometric mean of the smoothed and normalized 
spectra. The similarity of Si IV and C IV absorption profiles 
strongly suggests t > 1 and similar covering factor for both 
lines. Note the almost complete absence of C iv in emission. 
This is unlikely to be a result of absorption because the 
red-wing of the emission line appears smooth and without 
any features that can be attributed to absorption. The Lya 
emission line appears also weak but this might be attributed 
to an absorption of the red wing by N v and of the blue wing 
by Lya. Note also the prominent absorption features of the 
relatively low-ionization lines, in particular C ll and Al ill. 



* If no featureless window is found around ~ 1300 A, then fx at 
~ 1700 and 2100 A is used. If the 2100 A is unavailable (due to 
redshift), then two wavelengths in the region of 1700 A are used. 
^ The values used are from the VizieR catalog. Th ey do not ap- 
pear in the original paper of lSchneider et al.l fcoiOl) . 



Table 2. Parameters of the 8 objects that comprise the final 
BALQ data set. 



Name (SDSS J) 


z 




120954.14-1-142319.7 


2.7005 


-26.364 


145138.29-1-415401.0 


2.7425 


-26.908 


134934.14-1-245540.1 


2.7706 


-27.037 


143506.16-1-240144.8 


2.8096 


-27.717 


125224.35-1-144508.7 


2.9048 


-26.185 


111748.56-1-392746.2 


2.9125 


-27.578 


102321.90-1-493936.4 


3.0342 


-26.922 


134818.03+423205.1 


3.0661 


-26.938 



" The absolute magn itude in the SDSS i filter. Values are from 
ISchneider et al.l ll2010h . 

4.2 The non-BALQ control sample 

The procedure presented above requires a control non- 
BALQ spectrum to estimate the unabsorbed Lya proflle. 
We draw the control sam ple of non-BALQs from the SDSS 
DR7 quasar catalogue of (Schneider et all (|2010D . We select 
objects having values of z and Mi in the range of the BALQ 
data set (see above). We exclude objects th at are marked 
as BA LQ in at least one of the catalogues of'Scarin gi et al.l 
( 20091). I Allen et all l|201ll ') and Shcn ct al. (2011). This com- 
prises a sample of 1845 objects. Individual spectra are 
smoothed by a 22 pixel-wide moving- average filter and nor- 
malized using a mean f\ in the range of ~ 1600-1640 A (the 
same procedure as for the BALQ data set). The non-BALQ 
objects are sorted by the Lya emission strength. The mean 
normalized fx in a window of 1213.6— 1227.8 A, /a (Lya), is 
utilized as a proxy for the Ly a emission strength (the wave- 
length window is equivalent to —500 < v < 3000 km s~^ 
relative to Lya A1215.67). We choose the lower wavelength 
boundary to reduce the effect of intervening Lya absorp- 
tion. The upper boundary is chosen to reduce the contam- 
ination from the Nv emission line. Objects having mean 
/a (Lya) < 1 (i.e., Lya emission peak is lower than the con- 
tinuum at ~ 1620 A), or S/N in the i fllter smaller than 3 
are excluded from the sample. This criterion excludes 73 ob- 
jects, yielding a control sample that contains 1772 objects. 

The non-BALQ objects are divided into sub-samples 
based on the Lya strength. The first sub-sample contains 
the top 1 per cent of the objects (18 objects) in terms of the 
Lya strength, the second - 5 per cent, and the successive 
sub-samples contain 10, 20 percent, etc. of the objects. A 
geometric mean spectrum is calculated for each sub-sample 
and for the complete control sample. Figure 3] presents com- 
posite spectra of the 1 and 50 per cent sub-samples and of 
the complete sample. The composite spectra are overall sim- 
ilar, except in the strength of Lya and Civ emission lines. 
In particular, there is a similarity in the continuum level 
blue- ward of Ly a which appears not to be signiflcantly af- 
fected by intervening Lya absorption. Note that the BALQ 
composite sp ectrum is redder tha n the control composite 
spectra (e.g., iReichard et al.ll2003l ). 

4.3 Estimate of the metallicity 

The structure of this section is as follows. We begin by an 
outline of how different physical absorption settings can af- 
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Figure 4. Composite spectra of the 8 BALQ sample (tliick line) and control sub-sample composites: the top 1 per cent (thin line), the 
top 50 per cent (dashed line), and the whole control sample (i.e., 100 per cent; dash-dotted line). All spectra are normalized using a 
mean fx in the range of ~ 1600 — 1640 A. The laboratory wavelengths of the bluest line of each multiplet are indicated (the transition to 
the ground-state is indicated for Cll). Note the presence of prominent absorption features of low- ionization lines, in particular Cll and 
Aim; and the similarity between absorption profiles of Si IV and Civ lines. 



feet the measured column density (Sec. I4.3.1[l . This is fol- 
lowed by a description of the column density measurement 
procedure for four assumed physical settings (Sees. 14.3.21 - 
14.3.5(1 . Then, the procedures for estimating the normalized 
absorption profiles and the outflow metallicity are described 
(Secs. [473.6l and r4.3.7l respectively). Finally, the direct lower 
limits on absolute metallicity of the BALQ-sample objects 
are presented (Sec. 14. 3T^ . 

4-3.1 What sets the measured column density? 

The measured column density critically depends on the as- 
sumed physical setting of the absorber. If the profile is set 
by a foreground uniform absorber with a constant covering 
factor in velocity space, then the profile is governed by r 
only and the column density can be readily measured. If the 
profile is set by a velocity dependent CF, then r ^ 1 is pos- 
sible and the column density cannot be constrained. There 
can also be a combination of these two extreme settings. 
Thus, we have assumed several different physical settings 
of the absorbing gas and measured the resulting intrinsic 
r(Silv) for each of the 8 objects in the BALQ date set. The 
doublet nature of the Si iv line can be utilized to rule out the 
r ^ 1 setting for absorption profiles that have a deep dip 
in the blue part of the profile without a similar dip in the 
red part. The intrinsic r(Silv) is an output of the first step 
of the metallicity estimate procedure described in Sec. 13.31 
The different physical setting cases are described below. 

4-3.2 Case (i) - Non saturated absorber with a constant 
CF and a N(v) with a Gaussian profile 

The absorption profile is mostly set by 7V(ion) as a func- 
tion of V and we assume a Gaussian velocity distribution 
with a standard deviation b. Note that our goal is to fit 



the main bulk of the absorption rather than the exact pro- 
file with its small features. We follow the method to esti- 
mate the intrinsic r(Silv) that is described in lBaskin fc Laoil 
(2008), and is briefiy reviewed here. Each intrinsic syn- 
thetic absorption system has four free parameters: the ve- 
locity shift Wshift, Doppler broadening parameter b, global 
CF, and the ionic column density A*'(ion). The first three 
parameters are mostly set by the overall shape of the ab- 
sorption profile. The ionic column density is determined by 
the intensity of the absorption, but it has some degeneracy 
with the global CF. The intrinsic synthetic profile is trans- 
formed to an observed synthetic profile by convolving the 
former with the same smoothing function as the spectra. 
We do not convolve the intrinsic synthetic profile with the 
line-spread function (LSF) of the instrument, since the LSF 
width (~ 170 km s~^) is significantly smaller than the width 
of the absorption (> 3000 km s^^). The observed spectrum 
is corrected by the synthetic observed absorption profile that 
is produced by assuming a set of free parameters. The set 
is varied until the correction produces a reasonable recon- 
struction of a continuum and line emission. The success of 
the reconstruction is determined by eye. The Ly a profile is 
then corrected using the resulting b, CF and Af(Si^^), and 
the equality in Eq.[ll where the [iV(H°)/Ar(Si^+)]min is de- 
rived at the estimated A'^(Si^+). 

We prefer to use case (i) as the prime one. This is be- 
cause the fitting method used for this case, although not 
robust, is more certain than the other methods (described 
below). These methods require an a priori estimate of the 
unabsorbed spectrum in regions affected by emission lines, 
whose shape is not well constrained. We use cases (ii)-(iv) 
only for objects for which case (i) yields a lower limit on the 
metallicity that is higher than solar, since we want to check 
whether these cases also require Z > Zq. 
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4-3.3 Case (ii) - Non saturated absorber with a constant 
CF and a N{v) with an arbitrary profile 

The absorber is modelled by r(Silv;ii) and a global CF. 
The difference between this method and the previous one 
is that r(Silv;w) is calculated directly at each velocity bin, 
which we chose to be the instrumental width (~ 70 km s~^). 
We do not assume that N{Si^'^;v) has a Gaussian profile, 
but rather calculate it directly at each velocity bin assuming 
thermal broadening that is smaller than the bin width (i.e., 
6 ~ 10 km s^^). The column density at each velocity bin is 
calculated through the relation t{v) = ^^/i2Aion dN(v) /dv 
[hereafter, N{v) = dN{v)/dv x Av, where is the instru- 
ment bin width]. The main difficulty to execute this pro- 
cedure is the blending of the Siiv AA1394, 1403 doublet. 
One need s to de-bl cnd the contribution of the t wo lines at 
a given A. IJunkkarinen. Burbidge fc SmithI l|l983l ) suggested 
to overcome this difficulty by measuring the absorption pro- 
file beginning from its bluest part, where only the Si iv A1394 
contributes to the absorption. Then, using this measurement 
and the ratio of /12A of the two Silv lines, the expected 
Si IV A1403 absorption is calculated at a velocity that is 
shifted by the velocity difference between the two lines. If 
the observed absorption at this velocity is larger than the 
expected Silv A1403 absorption, an additional Silv A1394 
absorption is added at that velocity. The procedure is car- 
ried out throughou t the w hole absorption profile. Note that 
lJunkkarinen et al] (|l983f l did not include effects of CF in 
their method. These effects are accounted for in this study. 

The contribution of the blue and red line of the 
Si IV doublet is de-blended as follows. The calculation of 
7V(Si^+;u) begins from the bluest part of the absorption 
profile, assuming a contribution only from the Silv A1394 
line. We use the relation between the normalized absorp- 
tion profile / and r, I{v) = 1 - CF -|- CFexp[-r(u)] (the 
determination of I{v) is described in Sec. 14.3.61) . This rela- 
tion is utilized until the separation between Silv A1394 and 
SilV A1403 in velocity scale is reached (~ 1940 km s~^). 
Then, N(Q^'^ -jv) is calculated by the relation 

I{v) = 1 -CF + CFcxp[-r(u)] X exp[--r(u- 1940)/2]. 

The LyQ profile is calculated using the calculated 
r(Silv A1394;i;) and Eq. [H where [iV(H°)/iV(Si^+)]min is 
estimated assuming two cases. First, the whole outfiow is 
treated as a single photoionized slab, e.g. the outfiow is pho- 
toionized by photons entering it at a surface where v = Q. 
Second, each velocity bin is assumed to represent a sepa- 
rate photoionized slab e.g., a clumpy outfiow with clumps 
that do not cover the ionization source from each other (i.e., 
no self coverage). The latter case produces higher metallicity 
because [iV(H'')/iV(Si^+)]min is larger for lower N{^i^+) (see 
Fig- El - Intermediate cases produce intermediate metallicity. 
Note that the derived metallicity is almost insensitive to the 
precise choice of the velocity bin width Aii, because the un- 
derlying relation used in the procedure is between r(Lya) 
and r(Silv). 

4-3.4 Case (Hi) - A velocity dependent CF along the line 
of sight 

In this case, we assume that the absorber acceleration is 
along our line of sight to the emission source, and that the 



absorption profile is mainly set by a velocity dependent CF. 
We do not try to measure r(Silv) but rather estimate its 
lower limit. The measurement of CF(ii) and r(Silv;u) is 
conducted iteratively. First, we assume t{v) = 00 and cal- 
culated CF{v) through I{v) = 1 - CF(i;) + CF(«) exp[-r(t;)] 
beginning from the bluest part of the absorption profile. 
When the separation between Siiv A1394 and Siiv A1403 
is reached, we introduce the complementary part into the 
relation 

I{v) = {l-CF(?;) + CF(?;)exp[-r(u)]} 

X{1 - CF(?;') -f- CF(t;') exphr(t;')/2]}, 

where v' = v — 1940 km s~^. If no physical solution to CF(u) 
at a given v can be found (i.e., the predicted absorption by 
Siiv A1403 assuming r ^ 1 is larger than the observed one), 
we assume for this iteration step CF(t;) = (i.e., Siiv A1394 
does not contribute to the absorption at this v) and calcu- 
late t(v'), which is then used as an input for the next iter- 
ation step. The process is iterated until a physical solution 
to CF(i;) is found for all v, or when the maximum num- 
ber of iterations is reached. The Lya profile is estimated 
using CF(i;), r(Silv A1394; v) and Eq. [TH Here we assume 
the minimal value of [iV(H°)/Ar(Si^+)]mm « 2.5 (see Fig.[3l). 
Note that this case also checks whether the profile can be 
fitted by CF(ii) and r ^ 1. If it can, the algorithm converges 
after the first iteration. 

4.3.5 Case (iv) ~ A velocity dependent CF perpendicular 
to the line of sight 

In this case, we assume that the absorber acceleration is 
perpendicular to our line of sight. The procedure is similar 
to the previous case, except we use the following equation 
to deblend the profile: 

I{v) = l-CF(t;)-f CF(u)exp[-r(u)] 

-CF(u') + CF{v') exp[-r(u')/2]. 

This case also checks whether the profile can be fitted by 
CF(i;) and t > 1. 

4.3.6 Estimating the normalized absorption profile 

An estimate of a normalized Siiv absorption profile I (i.e., 
absorbed flux divided by the estimated unabsorbed flux) is 
calculated in order to apply procedures in cases (ii)-(iv) de- 
scribed above. We approximate the unabsorbed spectrum in 
the Siiv region by a constant value. The value is set by a 
mean f\ at the same range that is used to normalize the 
spectra (~ 1600 — 1640 A). This provides a good approxi- 
mation, as determined visually, for objects requiring the use 
of cases (ii)-(iv). This treatment under- predicts the absorp- 
tion in the red wing of the proflle, as it does not take into 
account absorption of the Siiv emission. Thus, the implied 
Lya absorption is also underpredicted. 

4.3.1 The metallicity estimate procedure 

The minimal required A^(Il''), implied by the measured 
iV(Si='+), is derived using CLOUDY models with solar 
metallicity (see Sec. 3). A curve of E versus U is calculated 
for the measured A;'(Si^+) (e.g.. Fig.©. The CLOUDY mod- 
els are interpolated using a 0.1 dex grid in E and U. The 
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curve is utilized to derive the allowable range of A'^(H'') val- 
ues. The minimal possible A'^(H'') is then found and used to 
reconstruct the observed spectrum in the Lyo? region. The 
photoionization models of a static slab should describe to 
a good approximation the photoionization structure of an 
outflowing gas. The observed outflow velocities of n < 0.03c 
do not alter considerably the ionizing flux when the latter 
is transformed to the gas rest-frame. 

The non-BALQ control sample is utilized to estimate 
whether the reconstructed spectrum in the Ly a region, us- 
ing solar metallicity, is consistent with a Ly a emission pro- 
file; and if it is not - to give a rough estimate of a lower limit 
on the metallicity. The excess EW of the reconstructed spec- 
trum is calculated relative to each control sub-sample in the 
region blue- ward of Lya. An excess EW > 10 A IS con- 
sidered as unacceptable. The value of 10 A is based on a 
rough estimate of the associated systematic uncertainties in 
this procedure. If the reconstructed spectrum is inconsistent 
with the top 1 per cent sub-sample, the minimal A''(H'') is 
scaled down by a given factor in the 1.5 — 100 range, and 
the fit procedure is repeated. The maximal factor which still 
produces an excess EW > 10 A is adopted as the object 
Z/Zq. 

4-3.8 Metallicity estimate results 

Four objects are consistent with solar metallicity, assuming 
case (i): SDSS J120954.14+142319.7, 145138.29+415401.0, 
143506.16-^240144.8 and 111748.56+392746.2. One object, 
SDSS J125224.35+144508.7, has an observed spectrum in 
the blue- wing region of Lya, where the flux density is al- 
ready larger than the top object of the control sample prior 
to making any correction. Thus, the object is excluded from 
the analysis. 

Three objects are inconsistent with solar metallicity as- 
suming case (i). We check whether the other cases allow 
lower metallicity than case (i). Cases (iii) and (iv) imply 
the same lowest lower limit on Z/Zq for these objects. We 
present for each object the case that yields the lowest ex- 
cess EW relative to the top 1 per cent control sub-sample 
for Z — Zq. Table [3] lists the parameters used in case (i) 
procedure, the resulting excess EW relative to the top 1 per 
cent control sub-sample, and an estimate of the lower limit 
of Z/Zq. The table also lists the Z/Zq for cases (iii) and 
(iv), the case that yields the lowest excess EW and its value. 

Figure [S] presents the results of cases (i) and (iii) mea- 
surement procedures for SDSS J134934.14+245540.1 . Case 
(i) is presented in panel (a). The Lya reconstructed spec- 
trum is compared to the top 1 per cent control sub-sample. 
The reconstruction produces a blue-shifted Lya emission 
profile which is inconsistent with the observations of non- 
BALQs. The measurement of the intrinsic absorption pro- 
file using the case (iii) procedure is presented in panel (b) of 
Fig. [SI and the reconstructed spectrum in the Ly a region is 
presented in panel (c). SDSS J134934.14+245540.1 is con- 
sistent with T(t)) S> 1 for all v i.e., its absorption profile can 
be completely fitted by CF(i'). Case (iii) implies metallicity 
consistent with solar for this object. 

Assuming t(v) ^ 1 for the remaining two objects pro- 
duces Si IV A 1403 absorption that is significantly stronger 
than observed. This requires to adopt a finite t{v) for 
cases (iii) and (iv). Figure [6] presents the results of the 



case (i) and case (iii) measurement procedures for SDSS 
J102321. 90+493936.4. The minimal CF{v) and the resuh- 
ing t{v) are presented in panel (c). The spectrum of SDSS 
J102321. 90+493936. 4 remains inconsistent with solar metal- 
licity for all cases. The lower limit for this object is Z > 
2Zq from case (iii) measurement. Figure [7] presents the re- 
suhs of cases (i) and (iv) for SDSS J134818.03+423205.1. 
Case (iv) applied to this object yields a metallicity that is 
consistent with solar. Thus, only the absorption is SDSS 
J102321. 90+493936.4 yields a robust resuh for higher than 
solar metallicity. 



5 DISCUSSION 

A mechanism which can produce outflows with a high 
metal/H abundance ratio is described in Sec. [2] Radiation 
pressure is mainly exerted on metals, which are coupled to 
protons and electrons by the Coulomb force. For a suffi- 
ciently large flux of photons, the acceleration of metals by 
radiation pressure is larger than the deceleration due to 
Coulomb friction, and the metals can decouple (i.e., run 
away) from the H gas, and produce a pure metal outflow. 
Several observed phenomena that may be explained by the 
proposed mechanism are describe below (Sec. I5.l"|) . We also 
discuss possible direct observational evidence for this mech- 
anism (Sec. 15. 2|) . and regions of the active nucleus that 
are suitable for the metal runaway mechanism (Sec. 15. 3p . 
Finally, several caveats to the proposed metal separation 
mechanism are described (Sec. 15.4]) . 



5.1 Extreme high velocity outflows 

The metal separation mechanism may be relevant to some 
observations which at face-value are inconsistent with the 
radiation-driven outflow scenario. It may explain the ex- 
tremely high outflow v elocities, v > 50, 000 km s ~ ^, fou nd by 



Jannuzi et al] (liggeP. iRodrfguez Hidalgo et al.l (|2007f ) and 
Rodriguez Hidalgo. Hamann fc HalJ ( 2011 ). where one of 



the objects, PC 0935+ 417, appears to have a relatively 
modest L/I/Edd ~ 0.2 (|Rodrfguez Hidalgo et al] l201ll )Pl 
The terminal velocity Wterm of a radiation-driven outflow is 
^ A4 X i/I/Edd^'csc, where M is the force multi- 
plier and Vcsc = V^'^Kcp is the escape velocity at the origin 



The ionic colu mn densities reported by 

iRodn'guez Hidalgo et aL I ll201ll ') for PG 0935+417 cannot 
be utilized to constrain the object metallicity by using the 
method described in Sec. [S] The authors provide only upper 
limits on 7V(ion) for the low ionization species [including N(}i^) 
that is estimated using the Ly a and Ly /3 absorption lines] and 
for JV(N*+). The Civ absorption cannot be used to constrain 
the metallicity, since the ob ject shows variable absorption 
llRodn'guez Hidalgo et al.l l201ll'l and C IV was not observed si- 
multaneously with Ly a and Ly /3. The remaining absorption line 
is Ovi, which is a high ionization specie line (requires ionizing 
photons with hu > 8.4 Ryd) and implies a very modest absolute 
lower limit of N{H°) = 0.02N{O^+) for logN{0^+) 15 - 17. 
The reported lower Umit on Af(0^+) > 8.2 X 10^^ cm~^ results 
in iV(H'^) > 0.2 X 10^^ cm~^ which is within the reported upper 
limit on N(H'^). Note that the measurement of all absorption 
lines except C IV is hindered by the intervening Ly a forest for 
this object. 
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Table 3. Parameters of the absorption system for the BALQ data set objects which are inconsistent with solar 
metallicity assuming case (i), and the lowest metallicity achieved assuming cases (ii)-(iv). 



Name (SDSS J) 


fshift 
(km s~-^) 


m 

(km s~^) 


CF 


7v(si3+ja 




eW3 

(A) 




CasJZl 


EWEl 

(A) 


M 


134934.14+245540.1 


-3800 


1000 


0.70 


15.34 


16.76 


49 


18 


(iii) 


7 


1 


102321.90+493936.4 


-3700 


1000 


0.90 


15.48 


16.83 


39 


5 


(iii) 


27 


2 


134818.03+423205.1 


-3000 


900 


0.80 


15.65 


16.87 


26 


4 


(iv) 


5 


1 



" Doppler broadening parameter. 

^ Log of the ionic column density in units of cm~^. 

Log of the minimal column density in units of cm~^. 

The excess EW relative to the top 1 per cent control sub-sample for case (i) and solar metallicity. 
^ Lower limit on the metallicity in units of Zq for case (i). 

f Both cases (iii) and (iv) imply the same lower limit on Z. The case that yields the lowest excess EW is indicated. 
9 The excess EW relative to the top 1 per cent control sub-sample for the indicated case and solar metallicity. 
^ Lower limit on the metallicity in units of Zq for the indicated case. 



of the outflow, where vk b t, is the Keplerian v elocity (e.g., 
IChelouche fc Netze dlioOll : [Laor fc Brandtll2002l ). The radia- 
tion pressure is exerted on the metals, which constitute only 
~ 1 per cent of the total mass for Z = Zq . A higher metal/H 
outflow implies the metals drag with them less H, and thus 
M oi: Z, and therefore Utcrm oc y/Z. Since Zq ~ 0.01, a 
pure metal outflow i.e., with Z ~ 1, will reach utcrm larger 
by a factor of 10. Only ~ 1 per cent of the total mass, which 
resides in metals, is being accelerated. Since the radiation 
pressure is unchanged, the acceleration is 100 times larger, 
and for the same acceleration region size, the final speed is 
10 times larger. 

The runaway me chanism may also explain the report 
bv lCottis et al.l (|201G| ) of the apparent lack of Ly q-N v line- 
locking signature in the Civ absorption trough in objects 
that meet all of the five physically motivated criteria to 
possess this signature, as suggested by Arav (1996). The 
high metal/H outfiow most likely can be approximated as a 
multi-component fiuid with a weak coupling between differ- 
ent components. The weak coupling is due to the relatively 
low number density of the fluid constituents. An ion runaway 
will cause a drop to ~ 10"'' in the particle number density, 
compared to a solar abundance gas. The coupling between 
an ion and an ion-fiuid due to Coulomb scattering is ~ Z^^-^ g 
(~ 4^ 20) times larger than that between an ion and a 
H-fiuid, where Ziou,a is the charge of the ion-fluid (note that 
■^ion fl = 1 in Eq. [HI since the equation is written for an am- 
bient fluid composed from protons and electrons only). This 
gives a ~ 1/(10~^ x 20) = 50 drop in the coupling force of 
an ion to an ion-fluid. Thus, the coupling between C and N 
is negligible in the outflow, and no Lya-Nv line- locking is 
expected for the C absorption lines. 

A direct observational manifestation of the weak cou- 
pling between different ions of a high metal/H outflow can 
be a difference in the absorption velocity distribution of dif- 
ferent ions, which cannot be simply explained by a change 
of ionization with the outflow velocity. Since, Wterm oc VM, 
M oc /i^(Aion)/i2/n^ion (Eq. m. Assuming the same illu- 
minating point source with fi, oc for all ions, yields 
i'term(ion) OC " /i 2 /"lion , aud oue needs to sum up all 
resonance transitions contributing to the acceleration of a 
given ion. It may be difficult to observationally identify rela- 
tively narrow extreme high velocity outflows, if different ions 



reach different utorm- The current technique to identify such 
systems is by looking for similar velocity shifts of the narrow 
absorption proflles, generally of C i v, N v and O vi, as for ex- 
ampl e measured in PC 0935+417 (Rodriguez Hidalgo et al-l 



5.2 The outflow metallicity 

We do not find robust evidence for a metal runaway (say 
by finding Z > IOZq) in the objects analysed here. The 
measured metallicities are similar to the a bsolute metal- 
licities found in earlie r studies of BALQs (IKo rista et ajj 
1 1996.) and NALQs (iHamann et all Il997l : lAraT et a l. 20o3; 
iHamann et al.|[201ll ). The moderate metal over abundance 
{Z > 2Zq) may refiect metal enrichment produced by en- 
hanced star formation rate close to the active nucleus. 

Utilizing relative metal abundances to test the runaway 
scenario is not straightforward. The relative metal abun- 
dances in the metal runaway scenario depend on both the 
initial metal enrichment, and the decoupling and acceler- 
ation of the different ions. Observed outflow s with over- 
abundance of some metals relative to other (e.g.. lWang et all 
[2009) may be caused by an SED which allows only the over- 
abundant metals to run away efficiently. In the local Uni- 
verse there is eviden ce for both sub- and su per-solar abun- 
dances in galaxies jZarritskv. Kennicutt fc Huchraj Il994l : 
ICallazzi et al.l |2005| ). observations at z > 1.5 point to- 
wards sub-solar abundances in central regions of ga laxies 
(|Cresci et al.ll201Gl : IJones et al.ll201Gl : lYuan et al.ll201ll ) , and 
the high metal/H outfiow may be characterized by a low 
metallicity metal abundance ratio. 



5.2.1 Prospects of future observations of nearby BALQs 

A high quality spectrum of the UV rest frame of nearby 
BALQs [z < 1) will allow to place a better constraint on 
the outflow absolute metallicity, as lower z AGN are signifl- 
cantly less absorbed by intervening Ly a absorption systems. 
The detected Lya absorption will reflect the intrinsic Lya 
absorption only, and will allow an accurate comparison of 
the relative Lya versus Siiv absorption profiles. Also, ex- 
treme high velocity metal systems will be detectable without 
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Figure 5. Spectrum reconstruction and absorption profile mea- 
surement for SDSS J134934. 14+245540.1 assuming cases (i) and 
(iii). Panel (a) presents the reconstructed spectrum assuming case 
(i). The laboratory wavelength of the blue line of hya and Si IV 
doublets is indicated by a vertical tick mark. If the reconstructed 
spectrum is assumed to represent the intrinsic unabsorbed spec- 
trum, then a very broad and blue-shifted Lya emission line is 
introduced. Panel (b) presents the measurement of the intrinsic 
absorption profile assuming a velocity-dependent CF. Procedure 
of case (iii) is carried out (see text) on the observed normal- 
ized profile (thick solid line), yielding an estimate of the intrinsic 
Si IV A1394 absorption profile (dashed line). The absorption pro- 
file of Si IV A1403 (dot-dashed line) is produced by shifting the 
Si IV A1394 profile by ~ 1940 km s~^. The reconstructed spectrum 
using the total Si IV absorption is presented (thin solid line). The 
unabsorbed emission is approximated by a constant (log fx = 0) . 
Note that the Si IV absorption can be fitted assuming t{v) 2> 1 
for all V. Panel (c) presents the reconstruction of hya. The ob- 
served spectrum (thick solid line) is corrected by the measured 
profile assuming solar metallicity, producing an estimate of the 
unabsorbed spectrum (thin solid line) which is compared to the 
top 1 per cent control sub-sample (dashed line). The two latter 
spectra are consistent i.e., the excess EW < 10 A. 

confusion by the Ly a forest, in particular since such systems 
tend to be narrow. 



5.3 Where can the high metal/H outflow 
originate from? 

There are two basic conditions for producing a high metal/H 
outflow. First, a metal runaway can be produced only if there 
is sufficient photon flux relative to the H density in a wave- 
length region that is absorbed by metal absorption lines e.g.. 
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Figure 6. Same as Fig. \E\ for SDSS J102321. 90-1-493936. 4 as- 
suming cases (i) and (iii). The intrinsic CF(v) (solid line) and 
exp[— t(d)] (dashed line) of Si IV A1394 are presented in panel 
(c). Note that the Si IV A1403 component predicts an absorption 
at D ~ 1000 km s"" (an emission feature in the reconstructed 
spectrum). This is expected because the normalized absorption 
spectrum is calculated disregarding the expected intrinsic Si IV 
emission. The Lyo reconstructed spectrum [panel (d)] is incon- 
sistent with the top 1 per cent control sub-sample (an excess EW 
of 27 A). Thus, indicating that the metallicity should be higher 
than solar. 



log/^/nn > -10.2 for T = lO'* K gas (see Eq. [T2j, which 
yields log [/ > 5 adopting the CLOUDY standard AGN 
SED. Second, the metal outflow can continue to accelerate 
only if the outflow does not become over-ionized after de- 
coupling from the H gas. This requires log U < 1 for the 
metal outflow, for the AGN SED, so that at least ~ 10"^ 
of C and N ca n absorb in Civ and Nv, respectively (e.g., 
lHamann|[l997l ). Consolidation of the two conditions requires 
a radiation-filtering shield that is further discussed below. 

The increase in U of the metal outflow is caused by two 
processes. 

(i) Once the metal outflow leaves its origin, the total num- 
ber density of free particles drops considerably. This leads 
to a large increase of U. The density drop is of the order of 
magnitude of ~ Z/20, where the factor of 20 represents the 
mean atomic weight of metals. Note that due to the lower 
number density the recombination rate is also lower (the loss 
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Figure 7. Same as Fig.|6]for SDSS J134818.03+423205.1 assum- 
ing cases (i) and (iv) . The resulting reconstructed Ly a spectrum 
assuming case (iv) is consistent with the top 1 per cent control 
sub-sample [panel (d)]; i.e., the spectrum is consistent with solar 
metallicity. 



of H means a loss of ~ lO'^ of the free electrons). The electron 
density drop leads to an increase in U by AU ~ 2 dex. 

(ii) The acceleration of the outflow gas causes an addi- 
tional decrease of the gas number density rigas- Since the 
thermal velocity of the outflow gas (6 ~ 10 km s~^) is 
much smaller than the outflow velocity {v ~ 10"^ km s~^), 
the expansion of the outflow perpendicular to its acceler- 
ation can be neglected, and the continuity equation yields 
rigas ccv^^ . Taking utorm to be ~ 30 x 10'^ km s~^, deflning 
the origin of the high metal/H outflow by w for which the 
metals run away {w ~ 10=* km s-^ Fig. [T]E1, and noting 
that the distance travelled by the gas until wtcrm is reached 
\r ~ |wLm/fflrad ~ (30 X 10*)VlO'' ~ 10^^ Cm; the value of 
Orad is adopted from Fig. [1] is probably much smaller than 
the distance from the ionizing source (r > 10^^ cm) i.e., 
the ionizing flux is constant throughout the outflow, yields 
a change in f/ oc nJaL oc i; of A?7 ~ logvtcrm/w — 1.5 dex. 



We use to rather than the typical velocity vtyp of the gas (e.g., 
a BLR cloud), since iityp is probably due to the Keplerian motion 
of the gas, and is perpendicular to the a^ad which is probably 
radial. 



A combined effect of the two processes is an increase of U 
by A(7 ~ 2 -f 1.5 = 3.5 dex. 

A shield, which filters the ionizing photons (i.e., low- 
ers U) without attenuating significantly the UV fiux that is 
absorbed by the metals, is needed to produce a metal run- 
away [e.g., the shield should produce log/^(UV)/nH > —10, 
T = 10 K and an effective log 17 < 0] ; and to prevent the 
resulting high metal/H outfiow to become over-ionized, and 
as a consequence to stall. The shield should fil ter the harder 
radiation that ca n over-ionize the meta l s (e.g..lMurrav et al. 
1995). Note that IChelouche fc Netze3 l|2003bl ') and lEveretl 
(2005) also invoke a shield while modelling a radiation- 
driven wind. They find that a shielded outfiow is accelerated 
more efficiently. 

There are two possible physical scenarios which can pro- 
duce a high metal/H outfiow without invoking an additional 
shielding component. 

(i) If the outfiow is launched from the BLR, then the BLR 
gas might provide the needed filtering. The outfiow might 
originate on the backside of the BLR cloud, behind most of 
the He""" ionization front for C^"^ and N'*'*' outfiows, and part 
of the H*^ ionization front for a Si'^"'" outfiow. The filtering is 
a result of the bound-free absorption edge, for example by 
He+ at 4 Ryd. The filter must allow production of N*+, 
and Si^^, and thus it should have r < 1 for photon energies 
E > 5.7, 3.5 and 2.5 Ryd, respectively. The filter also has 
to prevent further ionization of these species i.e., r ^ 1 for 
E > 7.2, 4.7 and 3.3 Ryd. Thus, a strong He""" edge, which is 
typically formed inside the BLR clouds, will form the right 
condition for the formation and acceleration of the C'^"'' ion. 

(ii) If the outfiow is launched from an atmosphere of an 
accretion disc (AD), then the accelerating radiation (UV) 
is found much closer to the gas, while the over-ionizing ra- 
diati on (EUV and above) is pr oduced at smaller disc radii 
(e.g-. IShakura fc Sunvae"^ll973i ). Both radiation sources are 
anisotropic. Thus, the gas can be illuminated by a suffi- 
ciently large to produce an ion runaway, while main- 
taining a small enough ionizing U that does not cause the 
outfiow to be over-ionized. For example, for a local black- 
body with temperature Tdisc, one requires that n„i that can 
over-ionize the ion (e.g., ionize C^^ to C'*"'') is smaller than 



for a Zp 







the specie number density i.e., < 10~ 
gas. We use ~ vf^{\t^)/huc, where Ath is the thresh- 
old wavelength required to over-ionize the specie. Inserting 
in Eq. [l2]for logT = 4 yields log /4Aion)//4Ath) > 9.1 
i.e., (i^ion/^'th)^exp[-/i(fion - t'th)/fcTdisc] > 10''-\ where 
Vi = c/Xi and the blackbody radiation is well approximated 
by the Wien law. This condition holds for Tdisc 35, 000, 
25,000 and 15,000 K for Nv, CiV and SilV, respectively. 
At these temperatures the local blackbody emission may be 
strong enough to produce and accelerate the ions, but too 
weak to destroy them. 

The metal outfiow over-ionization is readily prevented in 
the BLR scenario. The outfiow accelerates along an approxi- 
mately straight trajectory defined by the illuminating source 
and the BLR cloud, and the cloud continues to filter the radi- 
ation reaching the accelerating outfiow. If the metal outfiow 
is launched from an AD, it can be accelerated only along a 
trajectory that is approximately parallel to the AD surface. 
The contribution of the AD inner parts to the total radia- 
tion illuminating the outfiow increases as the outfiow reaches 
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larger altitudes above the disc surface. The radiation emit- 
ted by an AD becomes har der with a decreasing disc radius 
ijShakura fc SunvaevHToTi l. At a large enough altitude, the 
outflow is exposed to a hard enough total radiation that can 
over-ionize the gas, and the outflow stalls. 

5.4 Possible caveats to the metal decoupling 
scenario 

There are three possible caveats to the proposed metal de- 
coupling scenario that should be considered. 

(i) If the gas at the outflow origin is embedded in closed 
loops of magnetic field that wrap the gas, then the ions will 
move along the field lines unable to run away from the gas. 
However, if the magnetic field lines extend away from the 
gas, then the metal ions might be accelerated along them 
and run away from the gas, possibly producing an observable 
effect. 

(ii) The dynamical time for an ion runaway is ~ 
j-jth.o _ afj-iciw)] which is of the order of 5 min for 
the case presented in Fig. [T] [log /,^(Aion) = —4.2, lognn = 6 
and logT = 4]. If there are turbulent processes that mix the 
gas on time scales smaller than this, then the ion accelera- 
tion might be stalled, preventing the runaway. However, it 
is not clear what mechanism can produce nearly supersonic 
turbulence on the very small scales of 10^'' cm potentially 
involved with the acceleration. 

(iii) The effect of a bulk motion of the accelerated ions 
through the H gas is considered below. A bulk acceleration 
of ions should produce a current density that acts as an addi- 
tional deceleration forced However, if the accelerating ions 
are shielded by free electrons of the gas, a current density is 
not produced. The length scale of acceleration of a given 
metal ion can be approximated by /lacc ~ O.Svth.ion'^i^d- 
There is an electric shielding of the ion, if fea cc is much la rger 
than the Debye shielding distance, hu — kT/A-Knue'^, be- 
cause then the ionized fiuid (i.e., plasma) can effectively re- 
act to the ion motion. The condition of effective shielding 
is then 1 < ftacc/fto = 2.3 x W'' f-\f-''{Xi^,,)VT^ that 
can be approximated by log — 0.5 log T — 0.5 log nn < —4. 
This condition combined with the condition for a runaway 
(Eq. I12|) yields log tih — 3 log T < 4, which holds for a large 
parameter space (lognn < 12 and logT > 3). The shield- 
ing is effective while the ion velocity relative to the fiuid is 
w < Vth,c- For larger relative velocities {w > Wth.o) the bulk 
ions need to "carry" the shielding electrons. This results in 
a lower Coulomb frictional force at a given w, since the ef- 
fective charge of an ion and electrons it carries is smaller 
than Zion (see Eq. |8]). 



6 CONCLUSIONS 

We analyse the conditions required for a radiation-pressure- 
driven pure metal outflow, and possible observational evi- 
dence for its existence. We find the following; 

^■^ A current density produces a magnetic field. An interaction, 
between this magnetic field and charged particles that are coupled 
to the H gas, acts as a decelerating force on the bulk ions that 
produce the current density. 



(i) A metal ion embedded in gas subject to radiation pres- 
sure will run away when log fv (Aion) — log riH + log T > —6.2. 
(Eq. I12p . For an average AGN SED this can be converted 
to log U >5. 

(ii) To avoid overionization the gas must be protected by 
a shield, most likely in the form of a strong He"^ edge at 
4 Ryd. 

(iii) Photoionization models indicate for solar metallic- 
ity an absolute minimum of iV(H") = 2.5Ar(Si^+), implying 
r(Lya) ^ 1.8r(Silv A1394), regardless of U and E. 

(iv) Thus, a comparison of T(Si iv AA1394, 1403) and 
r(Lya), for different absorption geometries, can yield a di- 
rect constraint on the metals/H abundance r atio. 

(v) A search of the SDSS BALQ sample of lScaringi et al.l 
(|2009l ) yields a handful of possible candidates for an outfiow 
with a metallicity above solar. However, no robust evidence 
(i.e., Z > 10^0 ) is found for a metal runaway. 

(vi) High quality UV observation of lower z quasars, free 
of intervening Ly a absorption, can be used to obtain better 
constraints on the outflow metallicity. 

(vii) An ultra fast outflow of metals is expected to be 
produced if a runaway takes place. Since the metal ions are 
not expected to be collisionally coupled, different ions will 
likely have different terminal velocities. 

The physical conditions in AGN are deflnitely more 
complex than the simplistic scenario outlined in this paper. 
Clearly, models which simulate more accurately the vari- 
ous microphysics involved, including the caveats mentioned 
above, are required in order to get some quantitative pre- 
dictions on metal ions runaway in radiation pressure driven 
outflows in AGN. Such models can then be integrated into 
simulations of the large scale structure of AGN outflows. 
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